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I. INTRODUCTION 



In the past years increasing interest has been focused on semiconductor nanostructures.a This is mainly due to 
their low-dimensional character, which allows to tailor carrier quantum confinement as well as Coulomb interaction. 
As a result, this has allowed to fabricate nanostructured systems with a properly designed density-of-states whicbj 
in turn, exhibit an increased optical efficiency as well as a reduction of energy-relaxation and dephasing processes.S 
For the case of two- and one-dimensional nanostructures, i.e., quantum wells and wires, however, we deal with a 
partial carrier confinement, i.e., the single-particle energy spectrum is still continuous. This allows to describe their 
many-body ultrafast optical response in terms of the usual mean-field approaches, typical of bulk systems.H 

The real scientific and technological "revolution" in the field was the introduction of quasi zero-dimensional (OD) 
systems, called semiconductor quantum dots.u Compared to systems of higher dimensionality — like quantum wells 
and wires — they have a discrete, i.e., atomic-like, energy spectrum and, more important, they exhibit genuine few- 
carrier effects. Generally speaking, going from quantum wells and wires to quantum dots (QDsj) we move from 
many-electron systems to few-electron ones. This implies a radical change in the theoretical schemesci as well as in the 
experimental techniqueaJ used to study such quasi OD nanostructures, often referred to as semiconductor macroatoms. 
Apart from their relevance in terms of basic physics, these novel semiconductor nanostructures have attracted general 
attention because of their technological applications: these range from laserpemittera3 to charge-storage devices,S from 
fluorescent biological markersQ to quantum information processing devices.E3 

In QD's, the flexibility typical of semiconductors in controlling carrier densities has been brought to its extreme: it 
is possible to electrically inject-aiBgle electronaiil or to photogenerate in a QD a single Coulomb-correlated reieetron- 
hole pair, i.e., a single exciton.ll3ll3 It is even possible to detect the single-exciton decaying energy emission .113113 The 
quantized, atomic-like, energy structure of QD's allows for a rich optical spectrum and for a weak interaction of the 
QD system with environmental degrees of freedom (like phonons, plasmonSpfitc). This latter feature implies that the 
quantum evolution of the carrier subsystem is affected by low decoherenccEJ 

Moreover, their reduced spatial extension — up to few nanjometers — leads to an increase of two-body interactions 
among carriers and to stronger Coulomb-correlation effects.tl The latter may be used to design a variety of single- 
electron devices. In particular, as we shall show, they can be employed to design fully-optical quantum gates, as 
recently proposed in Ref. ^ Indeed, the continuous progress in QD fabrication and characterizationlla let us foresee 
a near future in which it will be possible to exactly tailor the few-carrier and optical properties of these OD systems. 
In this respect, a step forward has been recently made by the analysis and understanding of a single-QD excitoruc. 
emission spectrum,EjO that uncovered "hidden" symmetries in isolated QD structures, analogous to Hund's rulesE^I 
for real atoms. These symmetries imply that, under suitable conditions. Coulomb correlations among excitons in the 
same dot cancel. 

The primary goal of this paper is twofold. On the one hand, we shall present a detailed investigation of the electro- 
optical response of single- as well as coupled QD structures. More specifically, we shall focus on the combined effect 
of static electric fields and ultrafast multicolor laser pulses. Our investigation will present a variety of field-induced 
effects unexplored so far; in particular, we shall show how a properly tailored external field can be used to induce or 
reinforce exciton-exciton Coulomb coupling both in single and coupled QD structures. On the other hand, we shall 
discuss the application of such field-induced few-exciton effects to design a semiconductor-based fully-optical quantum 
information processing strategy.cB 

The paper is organized as follows: In Sect. || we shall introduce our theoretical approach for the analysis of the 
electro-optical response of QD structures; Section III presents a detailed investigation of the excitonic as well as 
biexcitonic response of prototypical semiconductor macroatoms and molecules in the presence of an applied static 
field; In Sect. IV our quantum information processing strategy is discussed and a few simulated experiments of basic 
quantum information/computation (QIC) operations are presented; Finally, in Sect, ^we shall summarize and draw 
some conclusions. 



II. THEORETICAL APPROACH 



The physical system under investigation is a gas of electron-hole pairs confined in a quasi-OD semiconductor struc- 
ture, i.e., a single as well as a multiple QD. In this case, the total Hamiltonian of our semiconductor nanostructure 
can be regarded as the sum of two terms, H = H° +H': A term H° describing the correlated electron-hole subsystem, 
i.e., free carriers plus confinement potential plus carrier-carrier Coulomb interaction, and a term H' describing the in- 
teraction of the carrier subsystem with coherent-light sources and environmental degrees of freedom, i.e., carrier- light 
plus carrier-phonon interactions. 
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A. Single-particle description 

Let us first consider the gas of non-interacting carriers, electrons (e) and holes (h) confined within the quasi-OD 

semiconductor structure. The quantum confinement can be described in terms of an effective potential Vc^^ whose 
height is dictated by the conduction/ valence band discontinuities. Since the energy region of interest is relatively close 
to the band gap Cgap of the semiconductors forming our heterostructure, we shall describe the bulk band structure in 

terms of the usual effective-mass approximation .0 In addition, since the confinement potential Vc^^ is jarfilowly- varying 
function on the scale of the lattice periodicity, we shall work within the "envelope-function" picture.t^ 

Within such approximation scheme, the non-interacting carriers in our quasi-OD structure are then described by 
the following Schrodinger equation 

A/jir) = e,/jVi/j(r) , (1) 

where m^/f^ is the bulk effective mass for electrons /holes while i/j denotes the set of single-particle quantum numbers, 
including charge as well as spin degrees of freedom.EHlHere, ipi/j^r) is the envelope function of state the eigenvalues 

e/h 

Ci/j correspond to the energy levels of the carriers induced by the confinement-potential profile Vc ; since the latter 
— for any realistic semiconductor nanostructure — is finite, the lowest part of the single-particle energy spectrum Si/j 
is discrete, while for increasing energies it evolves into a continuum. The different approaches commonly employed 
for the solution of Eq. (]l|) are described in App. according to the energy region of interest, they range from direct 
three-dimensional (3D) plane-wave expansions, to factorized-state solutions, or to simplified two-dimensional (2D) 
parabolic-potential models. 
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B. Coulomb-correlated carrier system 

Given the above single-particle representation {i} {{j}) for electrons (holes) — i.e., the set of 3D eigenfunctions 
'0i(r) = (r|i) (ipjir) = (r|j)) and the corresponding energy levels ei (ej) — let us now introduce the following creation 
and destruction operators for electrons and holes: 

\^)^cl\0)-.\0)^c,\^) , \j) ^ d]\0) \0) ^ d^\j) , (2) 

where |0) denotes the electron- hole vacuum state. Within such second-quantization picture, the single-particle Hamil- 
tonian, i.e., the Hamiltonian describing the non- interacting carriers within our OD confinement potential, can be 
written as: 

jjc = + H'' = ^ e4c, + E ■ (3) 

i 3 

The carriers (electrons and holes) within our quasi-OD nanostructure, however, interact via the two-body Coulomb 
potential V{y — r'). Due to such interaction, several correlation effects take place. Here, only processes conserving 
the total number of carriers are considered, thus Auger recombination and impact ionization are neglected. Such 
processes are known to become important only at very high densities and at energies high up in the bandO In this 
case the Hamiltonian describing carrier-carrier interaction within our single-particle z/j-picture can be written as 



2 / , ''»l»2«3*4 Si ^"12 '"is '"i4 



9 X/ ^'iJ 

32 ,j3 ,j4 

^ijij2i2 cijrfjjdjjCj^ , (4) 



2 "nnjsMdl^d^j^dj^dj^ 



n,i2 ji,j2 

where 



V,^,^i,i, = dr dv'rvS^Wvy)V{v - r')Vi2(r')V^;i(r) (5) 
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are the matrix elements of the Coulomb potential for the generic two-particle transition I1I2 l'il'2- The first two 
terms on the rhs of Eq. (jj) describe the repulsive electron-electron and hole-hole interactions while the third one 
describes the attractive interaction between electrons and holes. 

We stress the full 3D nature of the present approach based on the detailed knowledge of the 3D carrier wavefunctions 
Tp. The explicit evaluation of the above matrix elements for ageneric 3D confinement-potential profile Vc^^ , i.e., for 
a generic set of envelope functions ipi/j, is described in App. 

Combining the single-particle Hamiltonian in with the Coulomb-interaction term in we get the following 
many-body Schrodinger equation for our Coulomb-correlated system: 

H°|*) = (H^ + H'='=)|*) =£1*) . (6) 

Here, |^) denotes the interacting many-body state in our Fock space and £ the corresponding total energy. 
Let us now introduce the total-number operators for electrons and holes: 

Ne=5]cK, N, = ^4d^. . (7) 

i 3 

It is easy to show that the above global operators commute with the carrier Hamiltonian H° in (||) . We can therefore 
look for many-body states |^) corresponding to a given nuniber of electrons {Ne) and holes (Nh). In particular, we 
shall consider the case of intrinsic semiconductor materials,E3 for which N,, = N^; in this case the good quantum 
number is the total number of electron-hole pairs: N = Ng = Nh and the Schrodinger equation can be rewritten 
as: 



H°|Aa,) =£a„|Aw) , (8) 

where IAat) and denote, respectively, the A-th many-body state and energy level corresponding to N electron- hole 
pairs. 

For any given number N of electron-hole pairs we thus identify a subspace of the original Fock space, for which 
there exists a natural basis {I^at)}, given by the eigenstates of the single-particle Hamiltonian in (||): 

H^IM=e;«IM- (9) 

Here, In = Zi,i2, . . ■ ,«Ar; Ji,j2, ■ • ■ ,jN is a compact notation for the set of non-interacting electron and hole single- 
particle quantum numbers corresponding to our N electron-hole pairs. Indeed, we have: 

N 

\In) ^\{^njn})^ 114 A J^) (10) 

The non-interacting basis set in ( |lO|) constitutes the starting point of the direct-diagonalization approach used for 
the solution of the many-body Schrodinger equation (H). Indeed, we can expand the unknown many-body state IAat) 
over our single-particle basis: 

|Aa,)=^C/,^/|Z^). (11) 

In 

By inserting the above expansion into Eq. the latter is transformed into the foUowing eigenvalue problem: 

^^Hi,^-£xji^v^)u;:: = o, (12) 



where 



i?U = (^A.|H°|ZW)-Q„<5,„,^+Vvw (13) 

are the matrix elements of the carrier Hamiltonian H° in our single-particle basis. They are given by a diagonal — 
i.e., non- interacting — contribution plus a non-diagonal term given by the matrix elements of the Coulomb- interaction 
Hamiltonian in (^: V;^/^ — (Z^|Hcc|?^). Their explicit form — which involves the various two-body Coulomb matrix 
elements in (|^) — is given in App. ^ for the excitonic {N = 1) and biexcitonic {N = 2) case. 
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In the presence of Coulomb interaction, the Hamiltonian matrix in (^_3|) is non-diagonal; therefore, the interacting 
many-body states \Xn) are, in general, a linear superposition of all the single-particle states \In) [see Eq. (p])], 
whose coefficients U^^ can be regarded as elements of the unitary transformation connecting the single-particle to 

the interacting basis: U^'^ — {Ijs[\\n)- 

The numerical evaluation of our Coulomb-correlated states is thus performed by direct diagonalization of the 
Hamiltonian matrix H° in using a large — but finite — single-particle basis set. 



C. Interaction with coherent light sources 

The Coulomb-correlated carrier system described so far will interact strongly with electromagnetic fields in the 
optical range. For the case of a coherent light source — the one considered in this paper — the light- matter interaction 
Hamiltonian in our second-quantization picture can be written as: 

ti' = -mY.[^'^A4+^^^,d^C,] , (14) 

where E{t) is the classical light- field, and 

fJ'ij = fJ'buik / ^Ji{r)'tpj{r)dr (15) 



is the dipole matrix element for the ij transition, /ibuik being its bulk value. In the presence of a time-dependent 
coherent optical excitation the quantum-mechanical evolution of our electron-hole system will be described by the 
following time-dependent Schrodinger equation: 

^nJ^\^{t)) = H|*(t)) ^ (H° + H') \^it)) . (16) 

Contrary to the carrier Hamiltonian H° , the carrier-light term H' does not commute with the global number operators 
in (|^). Indeed, the two terms in (|lj) describe, respectively, the light-induced creation and destruction of an electron- 
hole pair. Therefore, N is no more a good quantum number and the many-body state at time t is, in general, a linear 
superposition of all the correlated A^-pair basis states: 

|*W) = EE«A«(^)|A^) . (17) 

By inserting the above expansion into the time-dependent Schrodinger equation ( p^ we get: 

iTi^axAt) ^ £x.ax,it) +J2J2^^^y .""^M) ' (18) 

where 

H'x^x'^, = (A^.|H'|AV) (19) 

arc the matrix elements of the light-matter Hamiltonian ( p^ ) within our interacting iV-pair basis {Ajv}- 

It can be easily shown (see App. ^) that the above matrix elements arc different from zero only for iV' = iV ± 1; 
this confirms that the only possible transitions are N^N + loiN + l^N which correspond, respectively, to the 
generation and destruction of Coulomb-correlated electron- hole pairs, i.e., excitons, discussed above. Moreover, we 
deal with well-precise spin selection rules: the only matrix elements in ( pj[ ) different from zero are those conserving 
the total spin of the carrier-light system. Indeed, the possible final states |A^) depend on the spin configuration of 
the initial many-body state |A^/) as well as on the polarization of the electromagnetic field E{t). In particular, we are 
allowed to create two excitons with opposite spin orientation (i.e., antiparallel-spin configuration) in the same orbital 
quantum state. In contrast, due to the Pauli exclusion principle, two excitons with the same spin orientation (i.e., 
parallel-spin configuration) cannot occupy the same orbital state. 

By treating Eq. (|l^ ) within the standard time-dependent perturbation-theory approach and assuming a monochro- 
matic light source of frequency w, we can define the absorption probability corresponding to the A^_j^ — > Aj^ transition: 

It describes the many-exciton optical response of our QD structure, i.e., the probability of creating a new exciton in 
the presence of — 1 Coulomb-correlated electron-hole pairs. 
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1. Excitonic absorption 



As a starting point, let us consider the so-called excitonic response, i.e., the optical response of our carrier system 
for the 0^1 transition. In this case, the initial {N = 0) state is the (electron-hole) vacuum state |0), while the final 
(TV = 1) state |Ai) corresponds to a Coulomb-correlated electron-hole pair, i.e., an exciton. Combining Eqs. (|T^) and 
(IH), for iV = 1 we have: 

|Ai) = Ef^/'^<4io)' (21) 

h 

where li = ii, ji denotes the single-particle electron- hole basis for iV = 1. 
The excitonic-absorption probability is then given by Eq. (po| ) with = 1: 

P^^{u;)^^\H',^,\'S{£x,-nu;) , (22) 

where 

F;^„-(Ai|H'|0) (23) 

is the matrix element of the light-matter Hamiltonian ( |l4|) for the — s- 1 optical transition. Its explicit form is given 
in App. 1^. 

The excitonic spectrum is finally obtained by summing the absorption probability in (|2^ ) over all possible final 
states |Ai): 

A-(..) = Ep,-(c.) . (24) 



2. Biexcitonic absorption 



Let us now come to the so-called biexcitonic response, i.e., the optical response corresponding to the 1^2 transition. 
In this case, the initial (A^ = 1) state coincides with the excitonic state |Ai) in (|l), while the final (A^ = 2) state IA2) 
corresponds to two Coulomb-correlated electron-hole pairs, i.e., a biexciton. Combining again Eqs. ( p^ ) and (p^), for 
A^ = 2 we get: 



(25) 



where I2 = iiji,i2j2 denotes the single-particle electron-hole basis for N — 2. 
The excitonic-absorption probability is then given by Eq. (EG) with N — 2: 



Tjbiex 



(^) 



27r 



(26) 



where 



H'x. = (A2|H'|Ai) 



(27) 



is the matrix element of the light-matter Hamiltonian ( p^ for the 1^2 optical transition. Its explicit form is given 
again in App. 

The biexcitonic spectrum is finally obtained by summing the absorption probability in (|2^) over all possible final 
states IA2): 



(28) 



We stress that, contrary to the excitonic spectrum in (p4[), the biexcitonic spectrum A''^'^^ is a function of the initial 
excitonic state X^. 



6 



Equations (|2^) and ( |26D will be employed in Sect. Ill to investigate the electro-optical response of single as well 
as coupled QD structures. However, for the case of ultrafast optical excitation and strong light-matter coupling, the 
above perturbation-theory picture can no longer be applied, and the time evolution of our many-body state \t/}{t)) can 
be obtained by solving the time-dependent Schrodinger equation in (p^. We stress that, contrary to the many-exciton 
absorption probability in (pO|), the number of excitons, i.e.. 

Nit) = (vI.(0|N,|vI/(t)) = (vI.(t)|N„|vI;(t)) , (29) 
is a continuous function of time and changes according to the specific ultrafast laser-pulse sequence considered. 



D. Interaction with environmental degrees of freedom 



Let us finally come to the interaction of the carrier subsystem with various environmental degrees of freedom, like 
phonons, plasmons, etc. They will not be treated explicitly; instead, we shall adopt a statistical description of the 
carrier subsystem in terms of its density-matrix operator 



the overbar denoting a suitable ensemble average. El Its time evolution can be schematically written as 



dt 



Pit) = ^p(t) 



H 



dt 



Pit) 



(30) 



(31) 



The first term describes the deterministic evolution induced by the carrier Hamiltonian H according to the well-known 
Liouville-von Neumann equation 



d , . 



H 



= ^[H,p(0] 



(32) 



while the second one describes a non-unitary evolutionj^l due to energy-relaxation and dephasing processes. The 
latter will be treated within the standard T\Ti model (see Sect. LVB). 



As for the case of the Schrodinger equation (p_6D, it is convenient to describe the density-matrix operator p — as 
well as its time evolution — within our Coulomb-correlated A^-pair basis. By combining Eqs. (O) and (BG), we get 



px,v(t)^a,(t)al,{t) : 



(33) 



the density matrix in the A-representation is bilinear in the state coefficients in (p^. 



E. The excitonic picture 



As discussed in Sect. II B, the generic Coulomb-correlated A-pair state lA^v) can be written as a linear combination 
[see Eq. (^l])] of the non-interacting electron- hole basis states in (p^). Such single-particle picture is used to compute 
Coulomb-correlated A^-pair states and energy levels via the exact-diagonalization approach described in App. |^. 
However, it is often convenient to adopt — in stead of a single-particle description — an excitonic- like picture, i.e., a 
quasi-particle number representation based on Coulomb-coupled electron-hole pairs. Aim of this section is (I) to show 
that, in general, such an excitonic description is not possible, and (II) to identify the basic requirements needed for 
such a quasi-particle number representation and therefore for QIC processing. 

To this end, let us introduce the following set of excitonic creation operators: 

|Ai)-a1J0), (34) 

where, as usual, |0) denotes the electron-hole vacuum state and Ai is the label for the generic excitonic (A^ — 1) state. 
By comparing Eq. ( |2l| ) with the above definition, we can write these excitonic operators in terms of the electron and 
hole operators, i.e., 
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Moreover, in view of the unitary character of the transformation [/, we get: 



44 = E^^'*<- (36) 

Ai 

If we now consider the exphcit form of the non-interacting basis states in (|o|), the generic A^-pair many-body state 
( |rT| ) can formally be written as 

I^^) = E^{'a"}I{^i}) (37) 

{Ai} 

with 

|{Ai}>=n4jO) . (38) 



The expansion in (37) would suggest to define a sort of excitonic number representation in terms of the iV-pair 
states |{Ai}). We stress that, in general, this is not possible. The point is that, in general, the set of states in ^3^ ) do 
not constitute a basis for our N-pair Hilbert subspace. This is intimately related to the fact that — contrary to electron 
and hole creation and destruction operators — the excitonic operators in Eq. ( ^5|) do not obey canonic commutation 
relations. In general, the commutator 



-Ai.A' 



^Ai 7 -''^1' 



(39) 



is itself an operator. This clearly prevents the introduction of number operators and, therefore, of a genuine quasi- 
particle number representation. 



As will be discussed in Sect. IV, two basic requirements are needed to perform quantum information processing: 
(i) the tensor-product structure of the "computational space" considered, and (ii) the SU (2) character of the rais- 
ing/lowering operators acting on our computational subsystems, known as "qubits" . The main question is thus to 
study if — and in which conditions — the coulomb-correlated electron-hole system discussed so far may act as quantum 
hardware, i.e., may be used to perform quantum information processing. This requires to identify a set of independent 
degrees of freedom, qubits, with a SU{2) character, the one of spin-^ systems. 

As a starting point, one should then check if there exist a set of independent excitonic degrees of freedom; this 
corresponds to verify that for any pair of excitonic states Ai and X'l the commutator of Eq. (^9|) is equal to zero. 

Let us now discuss the tensor product structure of our computational subspace. To this end let us consider again 
the case of two qubits a and b. Generally speaking, we know that the Hilbert space of a bipartite system is Ha d) Tib, 
where Ti-a/b &re the Hilbert spaces of the individual qubits. This means that if {|^a)} is an orthonormal basis set for 
Ha and {|4)} is an orthonormal basis set for Hb, then {\la) <E) \lb)} is a basis set for the whole computational space. 
What one needs to test is the possibility of writing the many-body ground state — corresponding in this case to a 
biexcitonic state A2 — as the product of two independent excitonic states A J and A^. This corresponds to verify that 

(A2|(|A?>|A?»=1. (40) 

Provided that the above requirements are fulfilled, let us now focus on the single qubit, i.e., on one of the independent 
excitonic states Ai. In this case, we want to check that the exciton creation/annihilation operators introduced in ( |3^ ) 
obey usual SU{2) commutation relations. More specifically, we are interested in defining the z-component pseudospin 
operator S^^ as 

^A, ^Ca„a, . (41) 

In order to check that this is really a z-component spin operator, we should verify that its average value over our 
many-body state is either plus or minus one. Deviations from this ideal scenario can be regarded as a measure of the 
leakage from our computational space due to the presence of external, i.e., non-computational, excitonic states. 



In Sect. [V we shall show that for prototypical GaAs-based quantum-dot molecules all the above requirements are 



well fulfilled and our excitonic system can indeed be used as quantum hardware. 
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III. ELECTRO-OPTICAL RESPONSE OF SEMICONDUCTOR QUANTUM DOTS 



In this section we shall analyze the electro-optical properties, i.e., the optical response in the presence of a static 
electric field F, of single as well as coupled QD structures. While the light-matter interaction is described by the 
Hamiltonian H' in (p^), the presence of a static field F can be accounted for by adding to the confinement potential 
Vc^^ in (|l|) the corresponding scalar-potential term: 

V;:'^{y) V^'^{v) ± eF • r . (42) 

Here, the ± sign refers, respectively, to electrons and holes. As discussed below, this sign difference will give rise to 
exciton-exciton coupling and significant field- induced energy renormalizations. 

Within the usual envelope-function picturejij the single-particle properties of our quasi-OD structure are described 
by the Schrodinger equation in (^. Similar to the case of semiconductor quantum wires,Eil a quantitative analysis of 
the whole single-particle spectrum requires a direct numerical solution of Eq. ; this can be performed using a fully 
3D plane- wave expansion, as described in Rcf. ^ and briefly recalled in App. |^. 

If, in contrast, our interest is limited to the low-energy range only, for most of the QD structures realized so 
far the carrier confinement can be described as the sum of two potential profiles, one acting along the growth (or 
perpendicular) direction and one affecting the in-plane (or parallel) coordinates only: 

y//''(r) = Vl'^ir^) + V'''^{v\\) . (43) 

As a consequence, the 3D carrier envelope function V'i/j can be factorized according to: 

V'./,(r)=V'i/,Jr^)V'!'„/,„(r||) (44) 

and the single-particle spectrum is the sum of the parallel and perpendicular ones: 

In this case, the original 3D problem is reduced to the solution of two independent Schrodinger equations, along 
the growth direction and within the parallel plane. This can be still performed employing the plane-wave-expansion 
approach described in App. ^ 

For most of the state-of-the-art QD structures we have strong confinement (few nanometers) along the growth 
direction, while the in-plane confinement potential V^^^ is much weaker. Moreover, as far as the low-energy region 
is concerned, the in-plane confinement is well described in terms of a 2D parabolic potential. For this case the 
Schrodinger equation within the 2D parallel subspace can be solved analytically — also in the presence of the static 
field F (see App. ^) — and thus the problem reduces to a numerical solution of the Schrodinger equation along the 
perpendicular direction. 

The analysis of the electro-optical response of semiconductor quantum dots presented in the reminder of this section 
is based on this parabolic-confinement model, whose derivation and validity limits are discussed in App. 

A. Single QD structure 

Let us start our analysis by considering the case of a single QD structure in the presence of an in-plane static 
electric field: F = (Fy, Fl = F, 0). Within the parabolic-confinement model previously introduced (see also App. 0), 

the quasi OD carrier confinement for both electrons and holes is described by an in-plane parabolic potential V^^'^ 



[see Eq. (A5)] plus a square-like potential V'^'^ corresponding to the interface band offset along the growth direction. 

As a starting point, we have considered an ideal QD structure characterized by the following material and con- 
finement parameters for electrons and holes: effective mass me = 0.05mo, mt = O.OSmo and parabolic-confinement 



energy fiuje = 40meV, huJh = 25meV; The well width is w = 50 A and the static dielectric constant [see Eq. (B3)] has 
been taken to be Eq = 12. Within this ideal QD model, the square- like potential profile along the growth direction 
is characterized by an infinite barrier height, i.e. Ae/h ^ 2m J^w'^ • choosing the above material and confinement 
parameters for electrons and holes, we deal with a very special case for which the set of electron and hole single-particle 



envelope functions i/j'^ and ■0'* coincide. Indeed, the in-plane spatial extension a in Eq. (All) is the same for electrons 
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and holes. Moreover, we shall discuss how this symmetry, not present in a realistic QD structure (see below), is 
related to special features in the optical response of the system ( "hidden symmetry" ) , as described in Ref . 

Due to the strong perpendicular carrier confinement, for both electrons and holes we deal with a single localized 
state; therefore, the low-energy single-particle spectrum is simp ly given by a sequence of equally spaced discrete levels 
corresponding to the 2D parabolic confinement [see Eq. ( |A12| )]. This scenario is not affected by the presence of the 
in-plane static field F, which manifests itself only through an overall red-shift A£ of the single-particle spectrum, 
known as Stark shift [sec Eq. ([A9|)]. 

In the absence of Coulomb interaction, both the excitonic and the biexcitonic absorption spectra [see Eqs. ( [2^ ) 
and (p^)] will exhibit optical transitions connecting the above single-particle energy levels. As usual, their amplitude 
is dictated by the corresponding optical matrix elements, i.e., oscillator strength, as well as by the combined state 
degeneracy, i.e., joint density-of-states (DOS). 

Figures H and g(B) show the excitonic-absorption spectrum for F = and F — 50kV/cm, respectively. Moreover, 
the excitonic-absorption spectra are compared to the single-particle ones, i.e., the ones evaluated in the absence of 
Coulomb interaction, which respectively correspond to the dashed curves in figure |] and to part (A) in figure |2[ As 
already pointed out, in this case the optical transitions connect the equally-spaced single-particle electron and hole 
states. 

As discussed in App. ^ for F = the only allowed optical transitions are those conserving the envelope function 
total angular momentum, i.e., m = —m' [see Eq. ( A20|)]; m oreover, due to the special symmetry between electrons and 



holes previously discussed, we have — n'^ [see Eq. ( |A1^ )]. Their amplitude is dictated by the joint state degeneracy, 
which for the single-particle case (see dashed line) is given by (n^ -t- 1). In contrast jfor finite values of the in-plane 
static field F (see part (A) in Fig. H), the above selection rules are relaxed (see App. ^ and we deal with new optical 
transitions corresponding to m + m' 7^ and ^ n^, not present in the field-free case. Moreover, in the presence of 
the static field the spectra exhibit a significant reduction in oscillator strength. This is ascribed to a reduction of the 



in-plane overlap between electrons and holes [see Eq. (A18)], due to the charge separation induced by the applied field 
[see Eq. (A8)]. This can be clearly seen in Fig. |[ where we show the single-particle electron and hole ground-state 



charge distributions (dashed curves) corresponding to the single-particle spectra of Fig. ||. For = (see Fig. ^) 
the electron and hole parabolic-potential minima coincide and, therefore, the two charge distributions exhibit the 
same symmetry center. In contrast, in the presence of the in-plane field F the two potential minima are shifted 
toward different directions. This, in turn, induces an electron-hole charge separation, as clearly shown in Fig. |^. Such 
charge displacement — which corresponds to the formation of an in-plane electrical dipole — is responsible for the 
oscillator-strength reduction in Fig. || previously discussed. 

Let us now come to the Coulomb-correlated case (see solid curves in Figs. |^, || and Figs. ^, U). In the presence of 
Coulomb interaction — which for the excitonic case (TV = 1) corresponds to electron- hole attraction — the main effect 
is a global red-shift of the Coulomb-correlated spectrum compared to the single-particle one. More precisely, for F = Q 
we find a relatively strong red-shift of the lowest optical transition (of about 20 meV) . For higher transitions this effect 
is reduced, which can be understood considering that high-energy states are characterized by an increasing spatial 
extension and, therefore, by a larger average distance between electrons and holes. Moreover, the Coulomb-correlated 
spectrum exhibits a transfer of oscillator strength toward low energies between quasi-degenerate optical transitions. 
This scenario is well established, and characterizes also systems of higher dimensionality, like quantum wells and 



For increasing values of the applied field we have a progressive reduction of the excitonic red-shift as well as 
of the oscillator-strength transfer, i.e., of the electron-hole attraction. This is confirmed by the inset of Fig. where 
the exciton binding energy is reported as a function of the applied field. 

In order to better understand the physical origin of this field-dependent behavior, we have carried on a detailed 
investigation of the excitonic wavefunction projected into the electron and hole subspaces. More precisely, by rewriting 
Eq. ( pi| ) in the coordinate representation, the two-body excitonic wave function is given by: 

vl/r(re,r^) = ^[/^V»(re)V',(r.) • (46) 

The square modulus of 5''^^ will then describe the conditional probability of finding the electron with coordinate re 
and the hole with coordinate v^. If we now integrate such quantity over one of the two coordinates we get 

/^(re) = 1 |vl/r(re,r,0|2dr, = Y.^^^^^' M^M^e) (47) 

ii' ,i 



and 



Hiyn) - 1 |vl'r(r„r,)|^dre = ^ C/^*t/,^.,V'*(r,0^,v(r,0 . 



(48) 
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The quantity f^/^ can be regarded as an effective single-particle probability distribution, which accounts for the 
electron- hole correlation described by the excitonic wave function in Eq. (^). In the absence of Coulomb correlation 

the transformation U reduces to the identity [U^^ = and the effective single-particle distributions fl^^ coincide 

with the square modulus of the single-particle wavefunctions of electrons and holes, i.e., /f (re) = |'(/'i(re)P, /^(r/i) — 

The effective charge distributions for electrons and holes -defined, respectively, in ( p7| ) and (Ey) — are plotted in 
Fig. ^ and in Fig. ^ for the ground-state-exciton case. As expected, in the presence of Coulomb correlation the 
charge distribution deviates from the corresponding Coulomb-free case (dashed curves). For F = (Fig. ^ the 
average distance between electrons and holes is very limited, which leads to a strong exciton binding (see Fig. |^). 
For increasing values of the applied field (Fig. ^ we see again an increasing charge separation. However, the effect 
is now reduced, compared to the Coulomb- free case (see dashed curves). This is due to the competition between the 



displaced parabolic potentials and the electron- hole Coulomb attraction (see also Fig. ^in Sect. IIIC). The latter is 
progressively reduced due to a significant increase of the electron- hole average distance (see again Fig. This also 
explains the reduction of the excitonic binding energy reported as inset in Fig. |^. 

The analysis presented so far suggests that the behavior of the system is governed by the following three charac- 
teristic lengths: 

(i) the radial extension a of the parabolic ground state, which in this case is the same for electrons and holes [see 



Eq. (All)] 



(ii) the excitonic Bohr radius 



a- = ^ , (49) 



where /j, is the reduced electron-hole mass; 
(iii) the total electron-hole displacement d = |d||' — djjl [see Eq. ([A^)]. 

Generally speaking, when a <^ a'^^ we are in the so-called strong-confinement limit: the carrier confinement is 
dictated by the single-particle parabolic potential only, which implies that the wavefunction of the excitonic ground 
state coincides with the product of the electron and hole single-particle wavefunctions, i.e., the expansion in Eq. ( |4^ ) 
contains just one term. In the opposite case, called weak-confinement limit [a 3> a"^^), the excitonic wavefunction 
depends on the relative coordinate only and resembles the 2D hydrogen-atom solution. 

For the case under investigation the situation is as follows. In the field-free case [d — Q), the excitonic Bohr radius 
{a^'' ~ 200 A) is of the same order of the electron and hole ground-state radial extension (a ~ 60 A) , which implies 
that the exciton wavefunction deviates from the product of the corresponding single-particle states. This is confirmed 
by the Coulomb-correlated carrier distribution of Fig. |3[ compared to the Coulomb-free one (dashed curve). However, 
we are not very far from the ideal strong-confinement limit previously discussed; Indeed, our numerical analysis has 
shown that the single-particle expansion in (^) can be limited to a relatively small number (6 x 6) of electron-hole 
states. For increasing values of the applied field — and, therefore, of the charge displacement d — the average distance 
between electrons and holes increases, thus reducing Coulomb-correlation effects. This is confirmed by the absorption 
spectra in Fig. |^ as well as by the carrier distributions in Fig. ^, where the difference between Coulomb-correlated and 
Coulomb-free results is significantly reduced. We can therefore conclude that the presence of an in-plane static field 
F induces a net electron-hole charge separation, which leads to a significant suppression of electron-hole Coulomb 
correlation. 



Let us now move to the biexcitonic response of our ideal QD system. As discussed in Sect. |II C| , contrary to the 
excitonic case investigated so far, the latter depends on the spin configuration of both initial {N = 1) and final {N = 2) 
Coulomb-correlated states. More precisely, due to the spin selection rules in the light-matter interaction Hamiltonian, 
we deal with two relevant cases only: the parallel- and the antiparallel-spin one. 

Let us consider first the antiparallel-spin configuration. In this case both excitons can occupy the low-energy 
orbital state. Figure ^ shows the biexcitonic spectrum (dashed curve) compared to the excitonic one (solid curve) 
for F — 30kV/cm. We can clearly identify a biexcitonic transition (see first peak of the dashed curve), which is 
blue-shifted with respect to the ground-state excitonic transition (see first peak of the solid curve). This energy 
renormalization is known as biexcitonic shifta 

AS = - £x, - £x' ■ (50) 
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This positive energy shift can be understood as follows: the applied field induces for both excitons the same charge 
separation (see Fig. ^ which results in a repulsive dipole-dipole coupling. This is confirmed by the field-dependent 
behavior of the biexcitonic splitting A£ reported as inset in Fig. ||. 

Moreover, as shown in the inset, in the field-free case, the dot behaves as an artificial atom and the energy ta 
add an exciton in a shell is, up to the first order in the Coulomb interaction, independent of the shell occupation.tll 
Indeed, within first-order perturbation theory, when F = the two excitons occupy the same spherically symmetric 
orbital ground state and for an ideal structure as the one we are considering, the biexcitonic splitting is exactly 
zero, because in this case the various attractive and repulsive Coulomb interactions cancel exactly.E3 This can be 
understood as follows. In the strong-confinement limit — which is not far from the regime considered here — and 
for antiparallel spins, the biexcitonic splitting can be very well approximated by the non- interacting single-particle 
probability distributions /q — |V'oP ^i^d Jq — IV'o P only, i.e., 

£q J J |r-r'| 

where A/ ~ /q ~ /o is the difference between the electron and hole single-particle probability densities. Due to 
the special symmetry of the ideal QD structure under investigation, in the field-free case we have A/ = and the 
biexcitonic splitting is zero as well. 

As already pointed out, in the QD structure under investigation we are not far from the strong confinement limit. 
However, since our calculation of the biexcitonic splitting is non-perturbative, we get, even in the field-free case, a 
non-vanishing AS. This small, but not negligible biexcitonic splitting (A£(F = 0) = 0.7meV, see inset in Fig. |^) 
measures the Coulomb-interaction contribution, underlying that the real ground-state biexcitonic wave function has 
contributions also from higher-level single-particle states. This value is compatible with the one given in Ref. |^. 
We shall now show that the above field-free behavior is due to the special choice of material and confinement parameters 
of the ideal QD structure investigated so far. To this aim, let us now move to the case of a realistic semiconductor 
macroatom. As prototypical system let us consider a GaAs/AlAs QD structure characterized by the following material 
parameters: effective masses rrie = 0.067too and ruh = 0.34mo, conduction- and valence-band offsets Ag = 1 eV and 
Ah = 0.58 eV, parabolic-confinement energies huje = 30meV and hujh = 24meV, well width w = 50 A, and static 
dielectric constant eq = 12.1. 

Figure ^ shows again the comparison between excitonic (solid curve) and biexcitonic spectrum (dashed curve) for 
the antiparallel case in the field-free case. As we can see, contrary to the result in Fig.j5| we now deal with a significant 
biexcitonic splitting AS also in the absence of the in-plane field (see inset in Fig. Indeed, for any realistic QD 
structure we deal with different spatial extensions cte and ah of the electron and hole single-particle in-plane ground 
states. In Fig. ^ we report the electron and hole single-particle charge distributions and (solid curves) as well 
as their difference A/ (dashed curve). As anticipated, due to the different material and confinement parameters, the 
charge distributions for electrons and holes do not coincide anymore. This, in turn, gives rise to local violations of 
charge neutrality, i.e.. A/ ^ 0, and therefore to a non- vanishing biexcitonic shift [see Eq. (|5l|)]. We finally stress that 
the presence of the in-plane static field leads to a further increase of A£ (see inset in Fig.W. 

Let us now move to the parallel-spin configuration. In this case we study the probability of creating a second exciton 
in the dot with the same spin orientation of the first one. Due to the Pauli exclusion principle, the two excitons are 
not allowed to occupy the same exciton state. As already pointed out [see Eq. (p8|)], the biexcitonic spectrum of 
the system depends on its initial excitonic state |Ai). In Fig. ||we compare the biexcitonic spectrum (dashed curve) 
with the corresponding excitonic spectrum (solid curve) for the field-free case. Here, the biexcitonic spectrum has 
been computed assuming, as initial state |Ai), the excitonic ground state. Let us focus on the low-energy part of 
the spectrum: as expected, due to the Pauli principle, the exciton ground state is forbidden to the second exciton, 
which can occupy any other high-energy state. Contrary to the antiparallel case (see inset in Fig. |^), we now deal 
with a negative biexcitonic shift A£, i.e., the lowest biexcitonic transition (solid curve) is red-shifted compared to the 
corresponding excitonic one (second peak of the dashed curve). As discussed in Ref. Ea, such energy renormalization 
(in this case of the order of 10 meV) can be ascribed to the attractive dipole-dipole interaction between the two excitons 
in the dot. Indeed, due to their different orbital quantum numbers Ai and X'l — and thus to their different spatial 
charge distributions — they exhibit significant Coulomb coupling. Contrary to the antiparallel-spin case previously 
discussed, (see inset in Fig. now the presence of an in-plane static field F leads to a reduction and eventually to an 
inversion of the biexcitonic shift (see inset in Fig. ^). This can be understood as follows: the application of the in-plane 
field leads to a progressive reduction of the attractive, i.e., spatially antiparallel, dipole-dipole Coulomb coupling; for 
high fields this is transformed into a repulsive, i.e., spatially parallel, dipole-dipole interaction and therefore to a 
positive biexcitonic splitting. The transition from red biexcitonic shifts to blue ones, occurs when the displacement 
induced by the electric field becomes of the same order or bigger than the excitonic Bohr radius. 
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B. Coupled QD structure 



Let us now consider the case of a semiconductor macromolecule, i.e., a coupled QD structure. In particular, as 
prototypical system we shall investigate the GaAs/AlAs coupled QD structure schematically depicted in Fig. ^. The 
material and confinement parameters are the same of the realistic single-QD structure previously investigated (see 
Figs. ^, ^, and |^): effective masses rric — 0.067mo and mh — 0.34mo, conduction- and valence-band offsets Ag = 1 eV 
and A/i = 0.58 eV, parabolic-confinement energies hwe = BOmeV and huih = 24meV, static dielectric constant 
So = 12.1. The square-like carrier confinement along the growth direction for electrons and holes is schematically 
depicted in Fig. ^ for our semiconductor macromolecule a + b. This is tailored in such a way to allow for an energy- 
selective creation/destruction of bound electron-hole pairs in dots a and b. Indeed, the width of wells a and b are 
slightly different, which corresponds to a blue-shift of about lOmeV of the single-particle optical transitions of dot b 
with respect to the corresponding transition in dot a. We stress that such energy shift is also present in the absence 
of inter-dot tunneling and Coulomb coupling. Moreover, the inter-dot barrier width {w ^ 50 A) is such to prevent 
single-particle tunneling and at the same time to allow for significant inter-dot Coulomb coupling. We stress that 
the geometrical and material parameters of the prqtpsed prototypical structure in Fig. ^ are fully compatible with 
current QD growth and characterization technologj^cJ. 

Let us discuss first the excitonic response of the semiconductor macromolecule (a + b) in Fig. |^. The excitonic 
(0 1) optical spectrum in the presence of an in-plane electric field F = 75kV/cm is shown in Fig. |l^. Here, the 
Coulomb-correlated result (B) is compared to the Coulomb- free one (A). The scenario is very similar to the single-dot 
case previously investigated (see Fig. |l|): for relatively strong values of the applied field, apart from a rigid red- 
shift, the Coulomb-correlated result is very similar to the Coulomb-free one. Here, the two lowest optical transitions 
correspond to the formation of direct ground-state excitons in dot a and &, respectively (see Fig. |^). In contrast, the 
high-energy peaks correspond to optical transitions involving excited states of the in-plane parabolic potential. Due 
to the strong in-plane carrier confinement — compared to the relatively large electron-hole charge displacement — the 
two low-energy excitonic states are expected to closely resemble the corresponding single-particle ones. 

Let us now come to the biexcitonic response of our semiconductor macromolecule. In view of the strong-confinement 
regime considered, we shall focus on the two ground-state excitons only. Moreover, since we are primarily interested 
in studying inter-dot Coulomb coupling, we shall consider the parallel-spin configuration. 

In Fig. |ri| the excitonic spectrum (solid curve) is compared to the biexcitonic one (dashed curve). The latter 
describes the generation of a second electron- hole pair in the presence of a previously created exciton (1 — > 2 optical 
transitions). In particular, here the previously generated exciton is assumed to be in dot a. As for the single-dot case 
previously investigated (see Fig. ^), the crucial feature in Fig. |ll|is the magnitude of the biexcitonic shift. For the QD 
structure under investigation we get energy splittings up to 8 meV (see inset in Fig. |ll|) . This can be ascribed again to 
the in-plane static field F, which induces, in both dots, the excitonic dipole previously investigated (see Fig. This, 
in turn, gives rise to significant inter-dot dipole-dipole coupling between adjacent excitonic states. The micioscopic 
nature of such exciton-exciton coupling is the same of the Forster process exploited by Quiroga and JohnsonEH for the 
generation of entangled states in coupled QDs. 

The physical origin of the biexcitonic shift AS is qualitatively described in Fig. ^ where we show the effective 
spatial charge distribution of the two electrons (ca and Cf,) and holes {ha and hb) corresponding to the biexcitonic 
ground state in Fig. |ll|. As we can see, the charge separation induced by the static field increases significantly the 
average distance between electrons and holes, thus decreasing their attractive interaction. On the other hand, the 
repulsive terms are basically field independent. This is the origin of the positive energy difference AS in Fig. |lj. 

Let us now inv estig ate the possibility of using such QD molecules as quantum hardware for QIC processing. As 
discussed in Sect. HE , to this end a few basic requirements should be fulfilled. First of all, the operators for the two 
ground-state excitons in dots a and b should commute. By evaluating the average value (over the biexcitonic ground 
state) of the commutator in Eq. (^9|), this came out to be negligibly small, thus confirming that these are indeed 
independent degrees of freedom. Moreover, due to the relatively large inter-dot distance — compared to the spatial 
extension of the carrier wavefunctions along the growth direction — the biexcitonic ground state in Fig. ^ is expected 
to closely resemble the product of the two excitonic states in dots a and b. Indeed, for the coupled QD structure 
under investigation we find that the scalar product in Eq. ( |40| ) gives a value of 0.99, very close to 1. The product 
structure for the bipartite system Hilbert space is therefore very well achieved. It is worth noticing that, in the case 
in which the two excitonic states are localized on the same dot, e.g., in the ground and first excited states, one gets 
a smaller value of about 0.9. This is a clear indication that the tensor product structure for the many-body state is 
much better achieved in a coupled QD structure than in a single QD system, as the one proposed in Ref. uB. 

Let us finally focus on the SU (2) character of our excitonic qubits. To this end, we have evaluated the average value 
(over the biexcitonic ground state) of the pseudospin operator introduced in Eq. (^) . By truncating the single- 
particle basis considering just the lowest energy level in each QD (strong-confinement limit), one gets (A2|5'|JA2) = 1, 
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thus confirming that the operators in Eqs. ( p9| ) and (|4j) are the generators of a SU{2) algebra. In contrast, far 
from the strong-confinement limit, we get a result which is of course dependent on how many single-particle states 
contribute to form an exciton. Therefore, if we calculate again the mean value of the commutator cons idering an 
enlarged single-particle basis set, we get deviations from the above ideal result. As anticipated in Sect. HE , this turns 
out to be a measure of the leakage of our qubit. 

We can therefore conclude that ground-state excitonic transitions in our coupled QD molecule fulfill all the basic 
requirements for a semiconductor-based implementation of QIC processing. They can be used as computational 
degrees of freedom, i.e., qubits, and the standard pseudospin language can be employed. 



C. A simplified model 



In this section we shall present a simplified model able to properly describe excitonic binding as well as interdot 
biexcitonic coupling. Its analytical solution will allow for an extremely quick way of identifying suitable parameter 
sets needed to employ the above coupled QD structure as semiconductor-based hardware for QIC processing (see 



Sect. IV) 



As a starting point, let us consider again the typical single-QD structures of Sect. Ill A, whose single-particle 
confinement is modeled in terms of a box-like potential of width a in the growth (or perpendicular) direction and a 
2D parabolic potential in the in-plane (or parallel) directions. As previously discussed, this allows a factorization of 
the original 3D single-particle problem into a perpendicular and a parallel one [see Eqs. ( ^4|) and (|4^)]. However, in 
the presence of Coulomb interaction such factorization is, in principle, no-longer possible. 

More specifically, let us consider the single-exciton problem {N = 1) in the presence of an in-plane electric field F; 
this can be described by the two-body Hamiltonian 

H = (r,) + H^r^) ^ . (52) 

As discussed in App . the single-particle Hamiltonians for electrons and holes can be written in the compact form 
[see Eqs. (^), (^, and (^)]: 

H^/' W = + Vl'\r^) + ^t\" + ^W'. • (53) 



Here, the presence of the applied field results in a displacement dn^'' [see Eq. ( A8)] of the parabolic-potential minimum 



as well as in a rigid energy shift ^tejh [see Eq. (A9)]. We want to show that for all the QD structures previously 
investigated Eq. (|5^) can be approximated to an analytically solvable form, and important quantities as wave functions 
or biexcitonic shifts can be easily estimated with a good degree of accuracy. 

In our QD structures the wavefunction is strongly confined along the growth direction by the square well potential, 
so that we can approximate (r\ — r'[)^ in the Coulomb term with its average value Z^. We choose I to be twice the 
average length related to the ground state of an infinite-height square well of width a, i.e., I = {2a / t:) [n"^ — 6)/12. 
It is thus possible to separate the Hamiltonian ( |52|) as H = Hy (re|| , r^jy) -I- H_L(re_L) -I- H_L(rftj_), where Hj_(rij_) — 
p1.^/2mi + V^{ri±) is the single-particle Hamiltonian along the growth direction — exactly solvable for the case of a 
parabolic potential as well as of an infinite-height square well. By further defining the center of mass (CM) and relative 
coordinates R = [me(rg|| -|- d|) + m/i(r/j|| — d|^)]/Af , (M = me + mh) and r — r/j|| — t^w , the in-plane Hamiltonian Hy 
becomes 

H||(R, r) - + \mujIR^ + ^ + i/i^'ld "'^ 



2M 2 2n 2' 



r ■ 



+ti{Lul^ u;l)Il-{d-r) 7^^, (54) 

where = mem?i/M is the reduced mass, = {l + A){u!'^+ujfJ/2,ujf = {1~ A){ll!^+ujI)/2, A — [{me — mh)/M]{uj'^ — 
'^D/i.^l + ^l) and 

d--d|j+d|f = eF(^^ + ^^) (55) 
denotes the total (electron-|-hole) field-induced in-plane displacement. 
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In the limit [ul — U}'j^)/{ui1 + u}\) ^ 1, the two coordinates are only weakly coupled, and the Schrodinger equation 
associated to the CM coordinate R is exactly solvable; in the general case, we shall concentrate on the ground 
state, though the generalization to higher states is straightforward. We can approximate the ground state of H|| 
as ^(r, R) = ^'j;(a;)x(j/, R), where x and y denote, respectively, the components of r parallel and perpendicular to 

the field F, x(y,R-) = e ^/(A^Tr)!/^ • e ^(A|j7r)i/^ = ^/WJ^r, and A^ = ^n/MujR. By averaging H|| over 
x(2/, R), we get the effective Hamiltonian "Reff — \h^r + ^'-^n. + + Veff{x), characterized by the effective 

potential 



Veff{x) = -^iui{x - df + Vc y-i^ j (56) 

with Vc{u) = —{e'^/e^/^^Xr)e'^Ko{u), Kq being the zero-order Bessel function. 

Since Ko{u) jlu ■ e^", in the limit x — s- cx), we regain the expected behavior for the Coulomb term 

2Xi e^/x^ + P e\x\ 

Notice that, considering the typical parameters of our systems {I ~ 20 A and Xr ~ 50 A), according to Eq. ( |57| ) 
there exists a relevant range of values for x where we cannot approximate Vc by its simpler asymptotic. Coulomb-like 
form. 

Since we are interested in the system ground state, we can approximate Veff around its minimum: 

Veffix) « + ^f^Co^x - xo)^ , (58) 

where Vq = Veffixo) and fiuj^ = d'^Vf.f f /dx^\^ . Within such approximation scheme, the eigenvalues and eigenfunc- 
tions of Heff can be evaluated analytically and, in particular, the approximate ground state eigenfunction becomes 

regime we are interested in (strong confinement and pronounced biexcitonic 
shift, i.e., large enough external field), the Coulomb attraction between electron and hole can be regarded as a per- 
turbation. In this regime, its main effect is to reduce the displacement d between electron and hole wave-function 
centers to xq, while the two single particle wave functions are, with a good approximation, rigidly translated. This 
can be understood by looking at Fig. |2|, where the potential Veff is plotted for three different values of the external 
field F. The solid lines correspond to the full potential, the dashed lines to its parabolic part, the dotted line to 
the Coulomb part (independent of F). For small and intermediate F the influence of the Coulomb field on the total 
potential is relevant. For intermediate fields the figure clearly shows that the minimum of the total potential is shifted 
with respect to the parabolic one. For strong and intermediate values of the applied field F, the effect of the shallow 
Coulomb potential on the region around the minimum of the total potential is mainly a rigid shift with respect to the 
unperturbed parabolic potential. For small fields, instead, the shape itself of the potential is definitely modified by 
the Coulomb term. 

In the regime of interest, we can write xo as 

xq = d~ Ax, (59) 

with Ax ^ d. By inserting ( |59| ) into dVef f / dx\xo = and considering, in the resulting equation, terms up to first 
order in Ax, the following analytical expression is obtained: 

Ax _ A, exp(0 Aif ^gP^ 



where ^ = (d^ P)/2Xl, Ki denotes the first-order Bessel function, AK = Ko{C) - ^i(C), A(A7<',i4:i) = 2llK + 
Ki{S^)/^, and a*^^ is the excitonic Bohr radius introduced in Eq. (^).. Notice that the prefactor Xr/a'^^ is a measure 
of the system confinement. In a similar way, setting uj = ujr + Auj in fj,oj — d'^Veff/dx'^\xo, we can calculate the effect 
of the Coulomb attraction on the potential shape: 



Alu Xr exp(^) /(i^ , , , , , 



uJr a"'-^ 20r VA2 
d A? 1 



2A{AK, Ki) ~ i [ak -f 2^^^^ 



3A(Ai^,ifi)n . (61) 
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In the strong-field limit \^/cP <^ 1 [see Eq. (p5|)], Aw/w^ = —Ax/d cx — (Ar/a'^^)(A^/d^), which shows that, in this 
regime, Coulomb corrections decrease very fast with increasing field. The condition Ax/d ^ 20% quantitatively 
defines the validity regime of the proposed approximation scheme. The latter coincides with the intermediate- and 
strong-field one, which is the regime of interest for the QD structures investigated below. It is also easy to show that 
in this regime the correction on the wave function due to Au / tOr is negligible with respect to the correction given by 



the shift Ax/d (see also Fig. 13) 



As previously discussed, the most important quantity for implementing our QIC scheme is the biexcitonic shift. 
This is in our c ase the energy shift due to the Coulomb interaction between two excitons sitting in neighboring dots 



(see Sect. [II B ). Within our model, we approximate the biexcitonic ground state as the product of two excitonic 
wavefunctions sitting in different dots and built according to the prescription given above. The wavefunction in 
the growth direction is approximated by a Gaussian of width 1/2 and the two dots are taken to have the same 
width a, i.e., the average of the two dots widths. This is reasonable since, for construction, the two dots have 
almost the same width. The desired biexcitonic shift A£ is then obtained averaging the corresponding two-exciton 
Hamiltonian over such factorized ground state. Within this approximation scheme, A£ becomes an easy-to-calculate 
sum of at most 2D integrals. In the corresponding validity region the model provides a good estimate for A£: Figure 
shows a comparison between the exact results (diamonds), the approximate result (solid curve) and the results 
obtained by neglecting completely Coulomb correlation in the wavefunctions (dotted line). The dashed curve shows 
the approximated results obtained setting Auj/ur = 0: as anticipated before, this correction is generally negligible. 

In order to implement our quantum computing scheme, system parameters as and F must satisfy some 

specific requirements. This determines the parameter space available in designing our QD structures. To this end, let 
us analyze the various constraints in details. 

First of all, (i) in order to have well-defined qubits, tunneling between dots must be suppressed; in agreement with 
state-of-the-art nanostructure technology, we have chosen for our calculations barrier heights of 1 eV (electrons) and 
0.58 eV (holes) and an interdot distance D = 100 A. On the other side, (ii) to implement our QIC scheme. Coulomb 
interaction between consecutive dots must be strong enough to produce a biexcitonic shift of the order of a few meV; 
this can be obtained either by tailoring in a suitable way the distance between the two dots or by varying the strength 
of the in-plane field F, since as a rough approximation 

Af (X , (62) 

where d is given by Eq. ( p5|) . Unfortunately, (iii) a side effect of a strong electric field F is to decrease the oscillator 
strength, and, accordingly, the system response to driving laser pulses; indeed, the electric field induces a spatial 
separation between electron and hole wavefunctions, thus decreasing their overlap (see Fig. If we now consider 
the confining parabolic potentials, (iv) in order to have well-defined quantum dots, the system must be in the strong- 
confinement regime previously introduced: the characteristic length associated to the parabolic potential in Eq. ( |54| ) 
must be smaller than the corresponding excitonic Bohr radius a^^ . On the other side, (v) as shown by Eq. ( |62[ ) and 
Eq. (^), a too strong parabolic confinement would in turn heavily decrease the biexcitonic shift A£. Last but not least, 
(vi) in order to be able to perform general QIC schemes, we must be able to energetically address specific excitations 
of the system unambiguously. This means that the peaks of interest in the optical spectra, namely ground-state 
excitonic and biexcitonic states, must be well-isolated from other high-energy transitions. This determines additional 
constraints on the value of fiuje and hujh ■ 

From the above discussion, it is clear that in order to satisfy at the same time all the requirements listed above 
[(i)-(vi)], the system parameters must be fine tuned so that a quick mean to scan the whole parameter space becomes 
necessary. The simplified model previously described came out to be quite efficient in performing such detailed 



analysis. The available parameter space for a reasonable field of F = 75kV/cm is shown in Fig. 14, Here, the typical 
error in the calculated values of A£ is 10-20%. The constraints imposed are A£ > 3.5 meV, oscillator strength greater 
than 0.15fibuik, T^^e > ^i-^/i, ^w/,. — AS > lOmeV, and Xr/a*^^ < 0.6. Based on this analysis, we have identified the 
parameter set used in the simulated experiments of QIC processing presented in the following section. 



IV. QUANTUM INFORMATION PROCESSING 

As anticipated in the introductory part of the paper, the advent of QIC§1 as an abstract concept, has stimulated 
a great deal of new thinking about how to design and realize quantum information processing devices. This goal is 
extremely challenging: generally speaking, one should be able to perform, on a system with a well-defined quantum 
state space (the computational space), precise quantum-state preparation, coherent quantum manipulations (gating) 
of arbitrary length, and state detection as well. It is well known that the major obstacle to implement this ideal 
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scheme is decoherence: the spoiHng of the unitary character of quantum evolution due to the uncontrollable coupling 
with environmental, i.e., non-computational, degrees of freedom. Mostly due to the need of low decoherence rates, the 
first proposals k)x experimental realizations of quantum information processing devices originated frcm specialties in 
atomic physics,L£l in quantum opticsJ23 and in nuclear and electron magnetic-resonance spectroscopy.!^ On the other 
hand, practically relevant quantum computations require a large number of quantum-hardware units (qubits), that are 
known to be hardly achievable in terms of such systems. In contrast, in spite of the serious difficulties related to the 
"fast" decoherence times, a solid-state implementation of QIC seems to be the only way to benefit synergically from 
the recent progress in ultrafast optoelectronics^ as well as in meso/nanostructure fabrication and characterization.Q 
Among the proposed solid-state ipiplementations one should mention those in superconducting-device physicsLa and 
in meso- and nanoscopic physics.Lj In particular, the first semiconductor-based proposal, by Loss and DiVincenzo, 
relies on spin dynamics in quantum dots; it exploits the low decoherence of spin degrees of freedom in comparison to 
the one of charge excitations. 

As originally eevisioned in Ref. ^ gating of charge excitations could be performed by exploiting present ultrafast 
laser technologyj23 that allows-to generate and manipulate electron-hole quantum states on a sub-picosecond time- 
scale: coherent- carrier- controlr3 In this respect, decoherence times on nano/microsecond-scales can be regarded as 
"long" ones. Based on this idea a few implementations have been recently put forward;L3 However, while in these 
proposals single-qubit operations are implemented by means of ultrafast optical spectroscopy, the control of two-qubit 
operations still involves the application of external fields and/or microcavity-mode couplings, whose switching times 
are much longer than decoherence times in semiconductors. It clearly follows that such proposals are currently out of 
reach in terms of state-of-the-art optoelectronics technology. 

As already pointed out in Ref. in order to take full advantage from modern ultrafast laser spectroscopy one should 
be able to design fully optical gating schemes able to perform single- and two-qubit operations on a sub-picosecond 
time-scale. Following this spirit, we have recently proposed the first all-optical implementation with semiconductor 
macromolecules .Ej 

Aim of this section is to review and discuss the semiconductor-based impl ement ation in Ref. whose quantum 
hardware consists of coupled QD structures, like those investigated in Sect. 



IIIB 



As described below, the crucial 



ingredient in our QIC scheme is the field-induced exciton-exciton coupling discussed in Sect. Ill, Indeed, the central 
idea in our QIC proposal is to exploit such few-exciton effects to design conditional operations. 



A. Quantum hardware and computational subspace 



As discussed in Sect. [IE, two basic requirements are needed for QIC processing: (i) the tensor-product structure 



of the quantum hardware, and (ii) the SU{2) character of the raising/lowering operators acting on the individual 



qubits. Based on the electro-optical-response analysis of Sect. [II, we can conclude that state-of-the-art coupled QD 
structures can be used as semiconductor-based hardware for quantum information processing. Indeed, as shown in the 
previous section, these requirements are well fulfilled by the prototypical QD molecules studied above. Our detailed 
investigation has shown that a proper tailoring of the QD confinement potential as well as of the inter-dot distance 
allows to identify a well- precis e subset of excitonic states, corresponding to intra-dot ground-state excitons. Indeed, 



as clearly shown in Sect. Ill B| (see Fig. ||), we can associate to each QD structure a ground-state exciton, i.e., its low- 



energy optical transition corresponding to the creation/destruction of a Coulomb-correlated electron- hole pair in that 
dot. We have shown that for these low-energy intra-dot optical transitions the corresponding exciton wavefunctions 
are localized in the various dots of the array; this allows us to label such subset of excitonic states acco rding to their 
host QD. In addition, in view of the relatively strong carrier confinement, leakage effects (see Sect. |IIE[ ) are expected 
to play a minor role. 

More specifically, following the second-quantization notation, let us denote with \n^) the absence (n^ = 0) — no 
conduction-band electrons — and the presence (n^ = 1) of a ground-state exciton — a Coulomb-correlated electron- 
hole pair — in dot z/; they constitute the single-qubit basis for the proposed QIC scheme: |0)^ and The whole 
computational state-space is then spanned by the basis set: 

\{n,}) ^(g,,\n,), (n, =0, 1). (63) 

The full many-body Hamiltonian H — H° -|-H' in ( p^ ) restricted to the above computational space will be described 
by the following matrix elements: 

HmK} - (M + H')| {<}) = ff{„„}{„;j + ^{„„}{<} • (64) 

They are the sum of two contributions: the first one is due to the Coulomb-correlated carrier-system Hamiltonian; 
the second is due to the carrier-light interaction Hamiltonian in (p^). As discussed in Sect. HC, the latter describes 



the creation/destruction of electron- hole pairs driven by ultrafast sequences of multicolor laser pulses. 
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Let us now focus on the term H° . As discussed in Sect. [IB, it preserves the total number of electron-hole pairs N ^ 
and this is still true within our reduced — i.e., computational — subspace. In general, the Hamilton Matrix H°^^_^ j^^, j 



is non-diagonal. However, for the case of the realistic coupled QD structure analyzed in Sect. [II B , non-diagonal terms 
are found to play a very minor role. In this case, the latter can be neglected and the Hamiltonian matrix H° is then 
diagonal in our number representation {n^}. This suggests to introduce corresponding number operators acting on 
our computational subspace: = \ni,)nt,{n^\ — The Hamiltonian H° reduced to our computational 

subspace can then be expressed in terms of such number operators. In particular, for an array of coupled QDs this 
can be written as: 

H° = ^ £^ + i ^ A£^^, n^, . (65) 



Here, S^, denotes the energy of the ground-state exciton in dot v while A£^,^i is the biexcitonic shift due to the 



Coulomb interaction between dots v and i/', introduced in Sect. HI [see Eq. ( pO[ ) and Fig. |ll[ . 

The effective Hamiltonian in ( |65| ) iias exactly the same structure of the one proposed by Lloyd in his pioneering 
paper on quantum cfiilular automata,EJ and it is the Model Hamiltonian currently used in most of the NMR quantum- 
computing schemes.EZI This fact is extremely important since it tells us that: 

(i) the present semiconductor-based implementation contains all relevant ingredients for the realization of basic 
QIC processing; 

(ii) it allows to establish a one-to-one correspondence between our semiconductor-based scheme and much more 
mature implementations, like NMR.eZI 

According to (|6^), the singlc-cxciton energy £1, is renormalized by the biexcitonic shift AS^,yi , induced by the presence 
of a second exciton in dot i/' : 

AE^yi Uyi . (66) 

In order to better illustrate this idea, let us focus again on the two-QD structure, i.e., two-qubit system, of Fig. ^ 
and fix our attention on one of the two dots, say dot h. The effective energy gap between |0)f, and |l)b depends now 
on the occupation of dot a. This elementary remark suggests to design properly tailored laser-pulse sequences to 
implement conditional logic gates between the two QD-qubits as well as single-qubit rotations. Indeed, by sending 
an ultrafast laser 7r-pulse with central energy hLO\j\no\ =£(,-!- AEhana-, the transition \n-b) — > |1 — rif,) (tt rotation) of 
the target qubit (dot h) is obtained if and only if the control qubit (dot a)|4s in the state \na). Notice that the above 
scheme corresponds to the so-called selective population transfer in NMRp3 alternative procedures used in that field 
can be adapted to the present proposal as well. 

Moreover, by denoting with IQ"^ the generic unitary transformation induced by the laser 7r-pulse of central frequency 
'^&['^a]j it is easy to check that the two-color pulse sequence U^U^ achieves the unconditional 7r-rotation of qubit h. 



B. A few simulated experiments 



In order to test the viability of the proposed quantum-computation strategy, we have performed a few simulated 
experiments of basic quantum information processing. To this aim, we have performed a direct time-dependent 
solution of the generalized Liouville-von Neumann equation in (|3^) restricted to our computational subspace, i.e., we 
have simulated the time evolution of the reduced density matrix: 

= • (67) 

As discussed in Sect. [I D| , this is governed by the total Hamiltonian 7i reduced to our computational subspace 
[see Eqs. (^) and (|65)] plus a non-unitary termed due to energy-relaxation and dephasing processes induced by 
envirouniental degrees of freedom, like phonons, plasmons, etc. The latter has been treated within the standard T\T2 
model.H 

We stress that the present density-matrix description, restricted to our computational subspace, does not ac- 
count for leakage effects, i.e., it neglects processes connecting states of the computational subspace to other — non- 
computational — excitonic states, and viceversa. Due to the strong-confinement character of our QD structures (see 
Sect. fU^ such leakage effects are expected to play a very minor role. A quantitative evaluation of the leakage dynamics 
would require the inclusion in our density-matrix description of non-computational states. 
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The above simulation scheme has been applied to the coupled-QD structure of Fig. ^ in thepresence of an in-plane 
static field F = 75kV/cm: Sa = 1.673 eV, Sb = 1.683 eV, and AS = 4meV (see inset in Fig. y). 

We shall start our time-dependent analysis by simulating a basic conditional two-qubit operation, the so-called 
controlled not (CNOT) gate. 

Our first simulated experiment is shown in Fig. The multi-color laser-pulse train (see central panel) is able to 
perform first a tt rotation of the qubit a: Then, the second pulse is tuned to the frequency St + thus performing a 
TT rotation of the qubit b since this corresponds to its renormalizcd transition energy [see Eq. (|66|)] when the neighbor 
qubit a is in state |l)a- The scenario described so far is confirmed by the time evolution of the exciton occupation 
numbers and nt (upper panel) as well as of the diagonal elements of the density matrix in our four-dimensional 
computational basis (lower panel). 

More specifically, at the beginning the system is in the state |0, 0) = |0)a (g) |0)b. Due to the first pulse at i = 
the computational state moves to the state |1,0) = |1)q (X) |0)fc. Finally, at time i = 1 ps the second pulse brings the 
system into the state |1, 1) = |l)a (S) |l)b. 

This realizes the first part of the well-known CNOT gate: the target qubit b is rotated if the control qubit a is 
in state \l)a- To complete it, one has to show that the state of the target qubit b remains unchanged if the control 
qubit a is in state |0)q. This has been checked by a second simulated experiment (not reported here) where the first 
pulse, being now off-resonant (with respect to dot a), does not change the computational state of the system. As a 
consequence, the second pulse is no more into resonance with the excitonic-transition energy of dot 6, since the latter 
is no more renormalizcd by the excitonic occupation of dot a. Therefore, the initial state of the system is |0,0) and 
the final one is again |0, 0). 

The simulated experiments discussed so far clearly show the potential realization of the CNOT gate, thus confirming 
the validity of the proposed semiconductor-based QIC strategy. However, the analysis presented so far deais,with 
factorized states, i.e., we have simulated the CNOT gate acting on basis states only. It is well knowncS that 

the key ingredient in any quantum-computation protocol is entanglement. Generally speaking, this corresponds to a 
non-trivial linear combination of our basis states. 

We shall now show that the CNOT gate previously discussed is able to transform a factorized state into a maximally 
entangled one. Figure ^ shows a simulated two-qubit operation driven again by a two-color laser-pulse sequence (see 
central panel). Initially, the system is in the state |0,0). The first laser pulse (at i = 0) is tailored in such a way 
to induce now a rotation of the qubit a: |0,0) (|0,0) + |l,0))/\/2. At time t — Ips a second pulse induces a 
conditional 7r-rotation of the qubit b: |0, 0) + |1, 0) |0, 0) + |1, 1). This last operation plays a central role in any QIC 
processing, since it transforms a factorized state ((|0)a + |l)a) ® lO)^) into an entangled state (|0)a ® |0)b -I- |l)a <8) |l)b)- 

As we can see, during the pulse energy-nonconserving (or off-resonant) transitionso take place; however, at the end 
of the pulse such effects vanish and the desired quantum state is reached. 

The experiments simulated above (see Figs. |l^ and 16) clearly show that the energy scale of the biexcitonic splitting 
AS in our QD nplecule (see Fig. ^l|) is compatible with the sub-picosecond operation time-scale of modern ultrafast 
laser technology .E3 



V. SUMMARY AND CONCLUSIONS 



We have presented a detailed analysis of the electro-optical response of single as well as coupled QD structures. 
More specifically, we have investigated the effect of a static electric field on the many-exciton optical response of quasi 
OD semiconductor nanostructures. Our analysis has shown that a proper tailoring of the single-particle confinement 
potential as well as of the inter-dot distance and applied field allows to induce and control intra- as well as interdot 
exciton-exciton coupling; this, in turn, may give rise to significant energy shifts of the optical transitions. 

This field-induced dipole-dipolc coupling constitutes the key ingredient of the proposed all optical implementation of 
QIC with a semiconductor-based quantum hardware. Our analysis has shown that energy-selected optical transitions 
in realistic state-of-the-art QD structures are good candidates for quantum-information encoding and manipulation. 
The sub-picosecond time-scale of ultrafast laser spectroscopy allows for a relatively large number of elementary 
operations within the exciton decoherence time. 

At this point a few comments are in order. First, we stress a very important feature of the proposed semiconductor- 
based implementation: as for NMR quantum computing, two-body interactions are always switched on (this should be 
compared to the schemes in which two-qubit gates are realized by turning on and off the coupling between subsystems, 
e.g., by means of slowly- varying fields and cavity-mode couplings); conditional as well as unconditional dynamics is 
realized by means of sequences of_ultrafast single-qubit operations whose length does not scale as a function of the 
total number of QDs in the array.Ej 
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Let us now come to the state measurement. In view of the few-exciton character of the proposed quantum hardware, 
the conventional measurement of the carrier subsystem by spectrahy-resolved luminescence needs to be replaced by 
more sensitive detection schemes. To this end, a viable strategy could be to apply to ouii-semiconductor-based 
structure the well-known recycling techniques commonly used in quantum-optics expcrimcnts.cHI Generally speaking, 
the idea is to properly combine quantum- and dielectric-confinement effects in order to obtain well-defined energy 
levels, on which design energy-selective photon-amplification schemes. An alternative approach would be to adopt a 



storage-qubit scheme, as recently proposed in Ref. 41 



The nanoscale range of the inter-dot coupling we employed for enabling conditional dynamics does not allow for 
space-selective optical addressing of individual qubits. For this reason, at least for our basic QD molecule (a + b), we 
resorted to an energy-selective addressing scheme. However, extending such strategy to the whole QD array would 
imply different values of the excitonic transition in each QD, i.e., S^, ^ £^/. This, besides obvious technological 
difficulties, would constitute a conceptual limitation of scalability towards massive Quantiim Computations. The 
problem can be avoided following a completely different strategy originally proposed by LloycEj and recently improved 
in Ref. by properly designed sequences of multicolor global pulses within a cellular-automaton scheme, local 
addressing is replaced by information-encoding transfer along our QD array. 

Finally, a present limitation of the proposed quantum hardware are the non-uniform structural and geometrical 
properties of the QDs in the array, which may give rise to energy broadenings larger than the biexcitonic shift. 
However, recent progress in QD fabrication — including the realization of QD structures in microcavities — will allow, 
we believe, to overcome this purely technological limitation. 
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APPENDIX A: EVALUATION OF SINGLE-PARTICLE PROPERTIES 

In this section we shall describe the numerical approach used for the evaluation of the single-particle properties — 3D 
wavefunctions and energy levels — for single as well as coupled QD structures. Within the standard envelope-function 
picturejl3 the non-interacting carriers in our quasi-OD structure in the presence of a static electric field are described 
by the Schrodinger equation (|l|) with the confinement potential in (|4^): 



- + V^"\v)±e-F -v 



^i/ji^) = e,;/jVj/j(r) • (Al) 



As for the case of semiconductor quantum wires,cJ a quantitative analysis of the whole single-particle spectrum 
ei/j requires a direct numerical solution of the above Schrodinger equation. This can be performed using a fully 
3D plane-wave expansion described in Ref. which is a straightforward generalization to QD structures of the 2D 
plane-wave expansion proposed in Ref. p4. 



As anticipated in Sect. Ill, when — as in this paper — we are interested in the low-energy range only, for most 



of the QD structures realized so far the carrier confinement can be described as the sum of two potential profiles 



acting along different directions [see Eq. ([43[)], which allows us to factorize the original 3D problem in (Al) into a 
perpendicular (_L) direction and a parallel (||) plane [see Eqs. (^J) and (|45|)]. Moreover, as far as the low-energy 
region is concerned, the in-plane or parallel confinement is well described by a 2D parabolic potential. In this case 
the Schrodinger equation within the 2D parallel subspace can be solved analytically (see below) and thus our problem 
reduces to a numerical solution of the Schrodinger equation along the perpendicular direction: 



2y72 



2^6/ /I 



vf\r^)±eF^r^ 



^t/sA^^) = ^t/,.^t/^A^^)- (A2) 



This has been solved using the plane-wave-expansion technique previously discussed. Within such approach, the 
unknown envelope function is written as a linear combination of plane waves, i.e.. 
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(A3) 



where G = are recip rocal lattice vectors associated to the periodicity box L. By substituting the above plane- wave 
expansion into Eq. (A2), the latter is transformed into the following eigenvalue problem: 

{Hgo' - e^Sac) bG' = 0, (A4) 



where Ti-cc £^re the matrix elements of the single-particle Hamiltonian in (A2) within our plane- wave basis. A direct 
diagonalization of TLgG' will then provide the desired perpendicular energy levels as well as the wavefunction 
coefficients ha- 

Let us now come back to the in-plane or parallel-subspace problem, which we treat within the 2D parabolic- 
confinement model previously mentioned, i.e.. 



The corresponding Schrodinger equation is of the form: 



1 



2m 



e/h 



K/hr\\ ± eF|| • r|| 



W I ■ (rii) = e / ■ w- I ■ (rii) 



(A5) 



(A6) 



It is well known that the presence of a static uniform electric field Fy does not change the parabolic nature of our 
confinement potential. Indeed, Eq. ( |A6| ) can be rewritten as: 



2me/h 2 ' " " 



V-" /■ (rii) = fe" - /^telh) ^Z-" /■ (i-ll) 



The presence of the applied field results in a shift 



= T- 



of the parabolic-potential minimum as well as in a rigid energy shift 



(A7) 



(A8) 



(A9) 



often referred to as Stark shift. We stress that in the presence of the electric field F|| we have different symmetry 
centers for electrons and holes; this, in turn, introduces significant modifications in the selection rules governing 
interband optical transitions (see below). 

As anticipated, the Schrodinger equation (A7) can be solved analytically. Due to the central symmetry of the 

problem (with respect to the parabolic-potential minimum d||^''), it is convenient to adopt a 2D polar-coordinate set. 

By denoting with r = |r|| — \ ^^'^ radial coordinate and with (p the azimuthal coordinate measured with respect 
to the field direction, we have: 

VL(r,(p) = a-(l"l+i), 

where 



TT (n + |m|)! 



(AlO) 



is the spatial extension of the harmonic-oscillator ground state — oj^/h - 

while >c!l"' {x) denotes the generalized Laguerre polynomial in the dimcnsionless variable x = r'^ / 



(All) 



being its oscillation frequency- 
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In view of the central symmetry of the problem, our quantum numbers are those of the angular momentum in two 
dimensions, i.e., a radial number n (n = 0, 1, 2, . . .) plus an orbital number m (m = —n, —n + 2, . . . , n — 2, n). 
The corresponding in-plane single-particle energy spectrum is given by: 

4lm = ^t^e/Zi (2n - |m| + 1) = huje/h [n^ + 1) , (Al2) 

where = 2n — |m| denotes the energy quantum number with degeneracy {n^ + !)■ 

The 3D single-particle energy spectrum is then given by the sum of equally-spaced energy-level sequences, i.e., 

= + eL» = + (^^^ + '^)^e/h ■ (A13) 



for each energy level e-^ — obtained by solving the eigenvalue problem in (A4) — we deal with an harmonic-oscillator 
spectrum with energy separation huj(./fi. 

Given the single-particle state factorization in (^), the corresponding dipole matrix elements in Eq. (^5|) can be 
factorized as well: 

M»j = A*b«ifc2'ij^2:|||j|| , (A14) 

with 

^ti^ = J ^t(^^)^tA^^)dr± (A15) 

and 

An= /<(r||)<(r||)dr„ . (A16) 
By inserting the plane- wave expansion (AS) into Eq. ( |A15| ), we get: 

^i..=E%^G ■ (A17) 



r,ip introduced in (AlO) 



Let us finally fo cus o n the in-plane integral in (A16). This can be rewritten in terms of the polar-coordinate set 



^L,n'm' = / ^jL{re,^e)^l^,{rh,^h)d^ . (A18) 



Here, r^/^ = jry — dy''''| and ^e/h ^re the corresponding azimuthal angles. In general, the two polar-coordinate sets 
for electrons and holes do not coincide. Indeed, in the presence of a static field Fn we have different symmetry centers 



see Eq. (A8)]. In contrast, for Fn — the two coordinate sets coincide (rii — r,(p ^ r^^ipe = rfnfh) and the above 



equation reduces to: 

^L.n'm' = I ^L(r, , rdrd^ • (A19) 



In this case — for which the symmetry centers for electrons and holes coincide — we deal with a number of well-known 
selection rules. In particular we have: 

-^nm.n'm' ^m+m' ■ (A20) 

This tells us that in the electron-hole generation process the total angular momentum is conserved. M oreov er, for the 



special case of equally extended electron and hole wavefunctions, i.e., ae = OLh (see discussion in Sect. [II A), we have 



'^nm^n'm' ^ ^ra+va' ^n^n' ^m+^ri' ^rie 5 (A21) 

i.e., the energy quantum number is conserved as well. 

In contrast, in the presence of the static field the above selection rules are violated, due to the fact that the 
Harmonic-oscillator wavefunctions in ( A10| ) are no longer eigenstates of the total angular momentum. As we shall 



see, the same considerations apply to the case of the two-body Coulomb matrix elements discussed in the following 
section. 
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APPENDIX B: EVALUATION OF TWO-BODY COULOMB MATRIX ELEMENTS 



In this section we shall describe the numerical approach used for the evaluation of the two-body Coulomb matrix 
elements. Starting again from the single-particle state factorization in (^^, the Coulomb matrix elements in Eq. (|^) 
can be rewritten as: 



(Bl) 



where 



,,^Jr|| -rj|) ^Jdr^l dr^^f (rl)F(r - r')^,t (B2) 

can be regarded as an in-plane effective potential obtained by integrating the original 3D Coulomb potential V 
— multiplied by the corresponding wavefunctions ip'^ — over the perpendicular direction. 

Let us consider the explicit form of the 3D Coulomb potential in (B2) written in terms of its Fourier transform 
along the perpendicular direction: 



e 

EaTT 



dqKo q 



So \r - r'l 

Here, Eq is the static dielectric constantjll q denotes the Fourier-transform parameter, while 

cos y dy 



Koix) = 



(B3) 



(B4) 



is the zero-order modified Bessel function. 



By inserting the Fourier expansion (B3) into Eq. (B2), we realize a factorization of the two space coordinates r± 



and r'j_. Indeed, by introducing the form factors 



the Coulomb matrix elements in (B2) can be simply written as: 



i; 1^211 



(ril 



e 



^0 (q 



(B5) 



(B6) 



Therefore, the evaluation of the effective Coulomb potential in (B2) reduces to the evaluation of the form factors T 
in Eq. (B5). To this end, by replacing the wavefunctions -0^ with their plane-wave expansion in (A3) we get: 

Tu'iq) - ^G^G'0(G' -G + q), (B7) 



GG' 



where 



Oik) 



3''"'^dr, 



(B8) 



are plane- wave overlap integrals over the periodicity region L, whose explicit form can be evaluated analytically. 

Therefore, for any shape of the perpendicular confinement potential, starting from the numerically computed 
eigenvectors be [see Eqs. ( [A^ ) and (A4)], we are able to obtain the various form factors T which, in turn, allow us to 
numerically evaluate the effective in-plane Coulomb potential in ( |B2| ) . Once the latter is known over a suitable space 
grid, the original six-dimensional integral in (^) is then reduced to the evaluation of the four-dimensional integral in 

(Bl). This requires some care, since the effective potential is singular for ry — rj| = 0. In order to eliminate such 
singularity, it is convenient to replace the integration coordinate rj| with the relative coordinate ry = r|| — r'| . Indeed, 
if we move to 2D polar-coordinate sets for the new integration variables ry and ry, the presence of the Jacobian 
function corresponding to the relative coordinate rn cancels the potential singularity. 
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We stress that — as for the case of optical matrix elements previously discussed (see App. ^) — in the absence of 
applied static fields the symmetry centers for electrons and holes coincide [see Eq. (A8)] and, due to global rotation 
symmetry, we get: 



(B9) 



and the numerical integration in Eq. (Bl) reduces to three variables only: on e an gular and two radial coordinates. 
As for the optical matrix elements previously discussed, the selection rule in ( |B9| ) describes the conservation of the 
total angular momentum in the Coulomb interaction process: mi + TO2 = m[ + mi, . 



In contrast, in the presence of an applied static field the selection rule (B9) is relaxed and we need to numerically 
solve the four-dimensional integral in Eq. (Bl). 



APPENDIX C: EVALUATION OF MANY-EXCITON STATES AND OPTICAL MATRIX ELEMENTS 



In this section we shall apply the exact-diagonalization approach introduced in Sect. II B to the excitonic {N — 
1) and biexcitonic {N = 2) case. Generally speaking, the method consists in a numerical diagonalization of the 
interacting-carrier Hamiltonian H° written in the single-particle basis {I^tv)} [see Eqs. (^, (pi]), and (p^)]. 

For the evaluation of excitonic states, i.e., states corresponding to a single Coulomb-correlated electron-hole pair, 
the proper basis set is given by the single-particle states in (noh with iV = 1, i.e.. 



|Zi)^|*iJi)=4^4J0) . (CI) 
The corresponding Hamiltonian matrix is given by: 

with 

= (*iJi|H1*'ij;) (C3) 

and 

h::,,,.^,^ = {iiJi\H''\Aj[} ■ (C4) 



Combining Eq. ( |Cl| ) with the explicit form of the non-interacting Hamiltonian H*^ in (^) and making use of the 
Fermionic commutation relations we get: 

In a similar way, combining Eq. ( |Cl| ) with the explicit form of the carrier-carrier Hamiltonian H'^'^ in (^, after a 
straightforward calculation we obtain: Hf^- = —Vi, /i' . 

Let us now come to the evaluation of biexcitonic states, i.e., states corresponding to two Coulomb-correlated 
electron-hole pairs. In this case the proper basis set is given by the single-particle states in ( p^ with N ^2, i.e., 

\l2)^\^lJl^2n)^4AAA^^) ■ (C6) 
The corresponding Hamiltonian matrix is given by: 

il'jli-ih-i'ij'Aj'-i " iljli2j2-i[j[i'2j'2 ~'~ iljli2j2;i[j[i2j2 \^ ' ) 

with 



and 



Hnn^.nA^i^^^^ = <^m^2J2[H^'KJ[^y^ ■ (C9) 
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Again, combining Eq. (C6) with the exphcit form of the non-interacting Hamiltonian in (Q) and making use of 
the Fermionic commutation relations, in this case we get: 

— {^ii + ^12 + £ji + £^2) ^iiji,i{j'2^iij2,i'2j'i 

— (fill + £12 + f-jl + £^2) '^iljl,«2-Jl'^*2j2,i'ij2 

+ (Cii + £^2 + ^jl + ^^2) ^iiji,i'2j'2^i2j2,i'ij'i ■ 

In a similar way, combining Eq. ( |C6| ) with the explicit form of the carrier-carrier Hamiltonian H"^*^ in 
straightforward calculation we obtain: 



H. 



1 

2 
+ 2 



{^ili2i[i'2 ~ ^i2ili[i'2 ~ ^n«2«2*'i ~'~ ^«2ili2*'i) {^hj2,j[j'2 ~ ^ilj2 j'aJl) 

^2 ~ ^2iljjj2 ^ ^'lj2i2il ^ ^j2jlj'2j[) ('^41 42 ~ '^^l 42 ^'1 ) 

^i2j2,i'2i'2 ~^ ^ilJli'li2^i2j2,i'2j[ "t" ^1 52 i'l '^i2il ,13 ~ ^1 J2 i'lij '^12^2 .i'l 
^i2j2,i{j'2 ~ ^iljli'2j2^i2j2,i[j[ ~ ^ilj2i'2j'i^i2jl,i[j2 ^" ^152*2 j2'^j2il,iiji 
^ilj2,i'2i'2 ^ ^i2jli[j2^ilj2,i'2j'i ~ ^i2j2i'^j[^iljl,i'20'2 ^^i2j2i'^j'2^iljl,i'20'i 



«1J1«2J 



^2^122^1 '^■'1^2, i'lia ^~ ^2il«2j2'^ili2,iii; + ^2 ^2 J'l '^'iJl .'i J2 ^2 52 «2:'2 " 



(CIO) 
after a 



(CU) 

Let us finally discuss the explicit form of the carrier-light matrix elements entering the many-exciton absorption 
probability in Eq. (20). 

For the excitonic absorption [see Eq. (|2|)] the corresponding matrix elements are defined in Eq. (p3|). Combining 
the explicit form of the carrier-light interaction Hamiltonian (|lj) with that of the generic excitonic state in Eq. ([2l|), 
we get: 



(C12) 



where U^^ is the unitary transformation from the non-interacting basis to the excitonic one, l\ is the non-interacting 
two-particle label corresponding to the single-particle states i\ and ji, while [ii^ = fii^j^ is the single- particle dipole 
matrix element given in Eq. (|l5|). 



For the biexcitonic absorption [see Eq. (^6|)] the corresponding matrix elements are defined in Eq. (p,T\). Combining 
the explicit form of the carrier-light interaction Hamiltonian ( p^ with that of the excitonic state in Eq. (^ij) as well 
as of the biexcitonic state in Eq. (Eq), we obtain: 



H'x2x. = -mJ2ui2 



{{K)2f^ln + K>Q i^Sm2.nn - 1) - ^;,,M:2,2 - ^;^2Mn..} ' 



(C13) 



where again is the unitary transformation from the non-interacting iV-particle basis to the interacting one, fi is 
the single-particle dipole matrix element, and I2 = iiji, «2j2 is the generic label for the non-interacting two-pair basis. 
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FIG. 1. Single-particle (dashed curve) and excitonic absorption spectra (solid curve) in the field-free case {F = 0). The 
inset shows how the exciton binding energy Af is reduced as the in-plane electric field F increases. 

FIG. 2. Single-particle (A) and excitonic absorption spectra (B) for an in-plane field F = 50 kV/cm. The excitonic spectrum 
in (B) — apart from a rigid shift due to Coulomb interaction — is now comparable to the single-particle one. Here, numbers 
from 1 to 4 identify corresponding transitions in each spectrum. 

FIG. 3. Effective electron and hole charge distribution for the ground-state exciton in the field-free case. The three curves 
correspond to non-interacting (n.i.) as well as to Coulomb-correlated e-h pairs as indicated. 



FIG. 4. Same as in Fig. ^ but for an in-plane field F = 50kV/cm. 

FIG. 5. Excitonic (solid curve) and biexcitonic optical response (dashed curve) for the antiparallel-spin configuration in the 
presence of an in-plane electric field F — 30kV/cm. The inset shows the biexcitonic splitting AS as a function of the in- plane 
field F. Notice that for F — 0, the latter becomes very small (A£^ = 0.7 meV), which is due to the special symmetry of the QD 
structure considered: Oe — Q-h (see text). 

FIG. 6. Excitonic (solid curve) and biexcitonic optical response (dashed curve) of a realistic QD structure for the an- 
tiparallel-spin configuration in the field-free case. The inset shows the biexcitonic splitting A£ as a function of the in-plane 
field F. Opposite to the symmetric case previously considered (see Fig. P), in this more realistic case the spatial extension 
for electrons and holes as well as their Coulomb matrix elements are considerably different; this is the physical origin of the 
positive biexcitonic shift in the field-free case (see text). 

FIG. 7. Effective electron and hole charge distributions for the ground-state exciton in the field-free case (solid curves) 
as well as their difference (dashed curve). Due to the realistic QD parameters considered, the charge neutrality is violated: 
A/ / (see text). 

FIG. 8. Excitonic (solid curve) and biexcitonic optical response (dashed curve) of a realistic QD structure for the parallel-spin 
configuration in the field-free case. The inset shows the biexcitonic splitting Af as a function of the in-plane field F. As can 
be seen from the spectra, the latter is now negative for F = Q. However, it becomes positive at high fields (see text). 

FIG. 9. Schematic representation of the electron and hole charge distribution as well as of the confinement potential profile 
in our Coulomb-coupled QD structure. The latter is tailored in such a way to allow for an energy-selective creation/destruction 
of bound electron-hole pairs in dots a and h. Moreover, the inter-dot barrier width [w ~ 50 A) is such to prevent single-particle 
tunneling and at the same time to allow for significant inter-dot Coulomb coupling (see text). 

FIG. 10. Excitonic response of the array unit cell (a+fe) in Fig. |^ for an in-plane field F — 75 kV/cm. The Coulomb-correlated 
result in (B) is compared to the Coulomb-free one in (A). 
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FIG. 11. Excitonic (solid curve) and biexcitonic spectrum (dashed curve) for an in-plane field F = 75kV/cm. Due to 
the well-defined polarization of our laser source, the structure in the biexcitonic spectrum (dashed curve) corresponds to the 
formation of an exciton in dot b given an exciton in dot a. One obtains a similar structure in the biexcitonic spectrum, 
symmetrically blue shifted with respect to the excitonic transition in dot a, if one considers as initial state an exciton in dot b. 
The biexcitonic shift A£ as a function of the in-plane field F is also reported in the inset. 

FIG. 12. Effective potential Veff{x) (solid line) as a function of the x coordinate for three different values of the external 
field F. Here, the following parameters have been used: me = 0.067mo, mn = 0.34mo, Tiuje ~ 30meV, huih = 24meV. The 
dashed line represents the parabolic part of 14// (a;) and the dashed-dotted line the Coulomb term. 

FIG. 13. Biexcitonic shift AS as a function of the in-plane field F. Here, the parameters used are the same as in Fig. |l^. 
The squares represent exact numerical results, the solid line the predictions of the model, the dashed line the predictions of 
the model after setting Au = 0, and the dotted line the results obtained by neglecting completely Coulomb interaction in the 
wave function. The inset reports the behavior of Acu/ur and Ax/d for the same range of applied fields. 

FIG. 14. Plot of the parameter space available for designing the QD molecule used as quantum hardware in our QIC 
implementation. This has been calculated using the proposed analytical model (see text). 

FIG. 15. Time-dependent simulation of a two-qubit operation realizing the first prescription (|1,0) |1, 1)) for a CNOT 
logic quantum gate on dots a and b (see text). Exciton populations Ua and Ub (upper panel) and diagonal density-matrix 
elements (lower panel) as a function of time. The laser-pulse sequence is also sketched (central panel). 

FIG. 16. Time-dependent simulation of a CNOT quantum gate transforming the factorized state |0,0) -I- |1,0) into a 
maximally entangled state |0,0) -I- |1, 1) for the coupled QD structure a -I- fe in Fig. ^ (see text). Exciton populations Ua and 
rib (upper panel) and diagonal density-matrix elements (lower panel) as a function of time. The laser-pulse sequence is also 
sketched (central panel). 
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